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Abstract. A coordinate system is constructed for a general accelerating observer in 
1+1 dimensions, and is used to determine the particle density of the massless Dirac 
vacuum for that observer. Equations are obtained for the spatial distribution and 
frequency distribution of massless fermions seen by this observer, in terms of the 
rapidity function of the observer's worldline. Examples that are considered include the 
uniformly accelerating observer as a limiting case, but do not involve particle horizons. 
Only the low frequency limit depends on the possible presence of particle horizons. 
The rest of the spectrum is 'almost thermal' whenever the observer's acceleration is 
'almost uniform'. 



1. Introduction 

Unruh [1] and Davies[2] seperately observed in th 1970's that an observer who 
accelerated uniformly through the Minkowski vacuum would see a "thermal" spectrum 
of scalar particles, whose perceived temperature depends on the acceleration. The 
same result holds for fermions [3, 4, 5, 6, 7, 8]. The constant acceleration case has 
been extensively investigated and reviewed [8, 9]. Most approaches to the constant 
acceleration case employ Rindler coordinates ([8, 10] are partial exceptions), relying on 
the symmetry properties of Rindler space, and often emphasizing the importance of 
the horizon. If the Killing symmetry of Rindler space or its associated horizon were 
indeed responsible for the result, then questions arise such as: would an observer who 
accelerates 'almost constantly' see a spectrum of particles which is 'almost thermal'? 
If an observer's acceleration is constant for a long but finite period of time, will that 
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observer see an 'almost thermal' spectrum during that period? If not, the stability of 
the Unruh result could be in question. 

Evidence from detector models [11, 12] supports the stability of the Unruh result 
for scalar fields. Sriramkumar [11] demonstrates that a detector which accelerates 
constantly, but which is turned on only for a finite time, will detect a spectrum of 
particles differing from the thermal result only by transients. 

Bogoliubov approaches to scalar particle creation have been considered in some 
classes of accelerating frames in 1+1 dimensions [13, 14, 15] and in 3+1 dimensions 
[16]. The method of Yang et. al. [15] was based on the 'instantaneous rest frame' of 
a given observer. Interesting results were obtained, but their interpretation is limited 
by the multivalued nature of the instantaneous rest frame [18, 20]. The 'instantaneous 
rest frame' will in general assign more than one time to any event further than l/g(t) 
from an accelerating observer - it does not foliate the spacetime. This leads in [15] to 
multivalued predictions for the particle content at any such event, raising problems of 
interpretation. The possible foliations considered in Sanchez [13, 14] and Whiting [16] 
are similar to those used here, and will be discussed further in the next section. 

Below we consider an arbitrarily moving observer in 1+1 dimensions, traveling 
through the massless fermionic Minkowski vacuum. The foliation we shall use is the 
"radar time" construction described in full elsewhere [17, 18, 19] and presented briefly 
in Section 2. In Section 3 this construction is used to derive an expression for the 
number of fermions measured by this observer in the inertial vacuum, and also the 
spatial distribution of those fermions. Examples are presented in Section 4, showing 
that observers who accelerate 'almost constantly' see an 'almost thermal' spectrum in 
the causal envelope of their period of acceleration, and further supporting the stability 
of the Unruh result. 

We use "natural" units throughout such that H — c — 1, and restrict our study to 
1 + 1 dimensions. 

2. Radar Coordinates 

The concept of radar coordinates for an accelerated observer is not new [18, 21]. They 
provide a set of coordinates having one-to-one correspondence with the entire region of 
the Lorentz frame with which the accelerated observer can exchange signals. Consider 
an observer whose worldline A in an inertial frame is described by: 



where T\ is the observer's proper time, and a(j\) is the observer's 'rapidity' at time T\. 
The exponential e a ^ is the obvious time-dependent generalization of the k of Bondi's 
/c-calculus [22, 23, 24]. The observer's acceleration is a(r A ) = da/dr\, and satisfies 
a 2 = —A^A^. The observer's worldline is completely specified by the choice of origin 
(i.e. 2^(0)) and the rapidity function a(r\), or by the choice of origin, the initial velocity, 
and the function a(r\). The radar time t(x m ) and radar distance p(x M ) are defined by: 




(1) 
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t = \(t + + t-) (2) 
P = -t~) 

where 

t + = time at which a signal from an event travelling in the observer's 
backward direction (i.e. t + z — constant) is received or sent. 

t~ = time at which a signal travelling in the observer's forward direction 
(i.e. t — z — constant) is received or sent. 

and are such that p is positive in the observer's forward direction and negative behind 
- this choice is convenient for 1 + 1 dimensions. From equations (1) and (2) it follows 
that t^ 1 are related to x ± by 

To see the connection between the r, p coordinates used here, and the t' , x' used in 
the scalar field work of Sanchez and Whiting [13, 14, 16], consider the most general 
transformation, written as equation (6.2) of [14] and described briefly in [16]. This is of 
the form: 

x + = g(x + ) and x~ = —f(—x ~) (3) 

where / and g are analytic functions and x r± = t' ±x' are the 'accelerated coordinates'. 
The coordinate t' in general has no connection with the proper time of any particular 
observer. In radar coordinates, the line p = is the worldline of the observer. The 
requirement that r coincide with the observer's proper time on this worldine gives: 

If £ = 1 <4) 

It is this condition that uniquely assigns the choice of foliation to the choice of 
observer, and which allows x^(t) to be written in terms of the observer's rapidity 
function, as in equation (1). Throughout [13], and in the main body of [14, 16], attention 
is restricted to the case / = g. Combined with the restriction (4), this leads to the time- 
reversal invariant case considered at the end of Section 3. 

The metric in radar coordinates is: 



ds 2 = e ("(- + )-°('-))(dr 2 - dp 2 ) 



The 'time-translation vector field' [17, 19] is simply k 11 -^ = J^, and the 
hypersurfaces S T are hypersurfaces of constant r. To illustrate these definitions we 
consider some examples. 
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2.1. Inertial Observer 

As a consistency check, consider an inertial observer with a velocity v relative to our 
original frame. Then a is constant, and x^(t\) = e ±a T\. The coordinates are now 
given by = e Ta x ± , so that 

1 / — a + . a —\ t VX y. , _ . X Vt 

t = |(e + e x ) = , p = |(e x^ — e x 1 



The radar coordinates of an inertial observer are just the coordinates of that observer's 
rest frame, as expected. 

2.2. Constant Acceleration 

The simplest nontrivial case is that of constant acceleration [20]. In this case a{r) = ar 
and we have x^(t\) = ±a~ 1 e ±aTA , which gives: 

T -Ya^i^l) ^MaV-a) (5) 

ds 2 = e 2ap (dr 2 - dp 2 ) (6) 



Region II 



Region U 



Figure 1. Hypersurfaces of simultaneity of a uniformly accelerating observer. 

These are Rindler coordinates, which cover only region U of Figure 1, as expected 
[18, 21]. The hypersurfaces of constant r are given by t TQ (x) = xtanh(aro). 

2.3. The 1 Smooth Turnaround 7 Observer 

Consider an observer (Alice say) with a(r\) = ar c tanh(r/r c ). Her trajectory is given 
by: 

d T e ±aTctanh(r/r c ) 

She has acceleration a(r) = acosh _2 (r/r c ), so is uniformly accelerating for \r\ « r c 
and is inertial for |r| >> r c . As the parameter r c — >• oo this case approaches the 
constant-acceleration case. However, there are no horizons for this 'smooth-turnaround' 
observer. 
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Figure 2 (A) . Hypersurfaces of simultane- 
ity for the 'smooth turnaround' observer 
in the (x,t) plane (in units of 1/a) for 
t c = 1/a. 



Figure 2(B). Hypersurfaces of simultane- 
ity for the 'smooth turnaround' observer, 
in the (x,t) plane (in units of 1/a) for 
t c = 3/a. 



Figure 2 shows the hypersurfaces of constant r for a smooth-turnaround observer (in 
units of 1/a) for r c = 1/a (A) and 3/a (B). In Figure 2 (A), in which the turn- around 
time is relatively short, we can distinguish four regions. In region F, where Alice is 
effectively inertial with velocity tanh(ar c ), her hypersurfaces of 'simultaneity' are those 
appropriate to her status as an inertial observer. The same applies in region P. In regions 
L, R however, her hypersurfaces of simultaneity are flat, and are contracted/expanded 
by the Doppler factors e aTc and e~ aTc respectively. This prevents the multi-valuedness of 
the 'instantaneous rest frame' construction [23, 24, 15]. (The point is discussed further 
in [18], in the context of the relativistic "twin-paradox".) Figure 2 (B) is dominated 
by the period of acceleration, and bears closer resemblance to the constant acceleration 
case of Figure 1. 



2.4- Acceleration at Late Times 
Consider an observer (Bob say) with 

a(r A ) =ar c log (l + e 7 ^) 
His acceleration is: 

n 



(7) 



(8) 



1 + e-W Tc 

Bob is stationary at early times, and his acceleration approaches a at late times. Bob's 
hypersurfaces of simultaneity are shown in Figure 3 (A) (in units of 1/a) for r c = 1/a. 
There is a future horizon at x~ = J °° dr (1 + e T / Tc )~ aTc , but no past horizon. 

Regions L, F are narrow in Figure 3 (A), and are bounded by the horizon. In 
the (p, r) plane, however, describing the space seen by the observer, the horizon is of 
course not visible (it corresponds to r~ — > oo for constant r + ), and the L,F regions 
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Figure 3(A). Hypersurfaces of simultaneity 
for an observer who accelerates at late times, 
in the (x,t) plane (in units of 1/a) for t c = 
1/a. 




Figure 3(B). Hypersurfaces of constant t 
(dotted) and constant x (hard) in the (p, r) 
plane (in units of 1 /a) of an observer who 
accelerates at late times, for r c = 1/a. 



are infinite in extent. This can be seen in Figure 3 (B), where surfaces of constant t 
(dotted) and constant x (hard) are plotted in the (p, r) plane. The values of t, x shown 
are {-10, -8, -6, -4, -2, 0, 2, 4, 6, 8, 60, 600} in units of 1/a. The line x = 600/a (top 
right) almost coincides with the line t = 600/a, each being almost null. The same 
applies to the lines x = 60 /a and t — 60/a. All of the (p, r) plane is covered by (x,t), 
although the lines of constant x (or of constant t) become exponentially tightly packed 
in regions F and R. The future horizon of Figure 3 (A) implies in Figure 3 (B) that the 
x = — A/a line (for instance) never crosses the t = —1/a line in the (p, r) plane. The 
point (x,t) = (— A/a, — 2/a) is outside the observer's causal envelope. 

3. Particle Density of the Dirac Vacuum 

The massless Dirac equation in flat 1 + 1 dimensional space can be written 



r 



dtp 
dx^ 







(9) 



and is frame-invariant provided that the 7 M -matrices transform as (contravariant) 
vectors. By choosing the flat-space Dirac matrices 7 = ( ^ J \ and 7 1 = ( ^ ^ 

(which satisfy {7^,7^} = 2n^ v ) and introducing the basis spinors 0+ = 0_ = (°), 
we can write the general solution as 

i) = i) + (x-)<p + + i,_(x + )<p_ (10) 

To define the (inertial) Minkowski vacuum, introduce the plane wave basis: 



0+, 



u p ,b,± = e 



^fipx^ 



<f>- 



for p > 



(11) 
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These modes are normalized to 2n5(p' — p), and have been decoupled into "forward- 
moving" modes (denoted with a subscript "F"), and "backward-moving" modes (with 
subscript "B"). 

The Dirac field operator ifj(x) can now be written: 

fa) = [ 7^ J2 i u p^+ a p^ + u p^- b U ( 12 ) 

^° <j=F,B 

and the Dirac vacuum |0m) can be defined by the condition a p>a \0M) = = &p, ct |0m) for 
all p, a. 

Consider now an arbitrarily moving observer, with coordinates (r, p) . With 
7 ° = |L7 + g^ 1 and 7 1 = §f7° + g^ 1 , equation (9) becomes 

Since 7°7 1 = 5' 00 7°7 1 , this can be written as 

^ = -nV|^Arf (13) 

Following [19, 25] we construct particle/antiparticle modes for this observer, by 
diagonalising the 'second quantized' Hamiltonian: 



H(t) 



where M 



- a ( T ^(I + 7°7 1 ) + e a(r+) (I-7°7 1 ) 



(14) 

(15) 
(16) 



We now seek the spectrum of the 'first quantized' Hamiltonian H\(t) = \(H n h + 
H^), which acts on states defined on E r . The inner product expressed on such 
hypersurfaces can be written in radar coordinates as: 



mhr = J W^ f 7°7V= J dp ^M<\> 

— CO — oo 

To deduce H^ h , consider 

oo 

(iW|0>E ro = \ J dp (H nh ^y M<\> 

— oo 

oo 

e) f { dJ^f 

= (^-za—l^ - J dp^M (M-ha— 



(17) 



;is) 
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where we have discarded surface terms at p — > ±00, and have assumed that the functions 
U + and U~ are continuous. This gives our Hermitian Hamiltonian Hi as 

#i(ro)= ~i*§- p 

+ % - [a(r + p)(l- ft) - a(r - p)(l + ft)} (19) 

The positive/negative eigenstates of this Hamiltonian will determine the parti- 
cle/antiparticle content of the Minkowski vacuum measured by our accelerated observer. 
By writing these eigenstates as ip T0 = + /-</>-, equation (19) gives 

uf+ = - i-r^ - 2 a ( r o - P)f+ 
df 1 

ouf-=i-±- + -a(T + p)f- (20) 
and since (jl^j = {-^+) T = — (a7^) r > these equations can be rewritten as 

t d f+ 1 ( -\f 
0JJ+ — 1- a(T ) f+ 

J+ dr- 2 v !J + 

= ^ + \^ + )f- ( 21 ) 

with solutions 

/ + = A^ T ^l 2 e~ WT ~ 

f_ = J B w e- Q(r+)/2 e-^ r+ (22) 

where A u and B u are normalisation constants. Comparison with equation (10) (and 
observing that r + is a function of function of x ) reveals that in (21) we 

have found not only eigenstates on E T , but have derived solutions to the Dirac equation 
(9) which are eigenstates for all r. This is in general only possible for the massless 
case, and is a consequence of conformal invariance. By choosing {A — 1, B — 0} and 
{A = 0, B = 1} respectively, we can define "forward-moving" and "backward- moving" 
modes: 

i/>u,,f,± = e a(T ~ )/2 e^ r >+, VW?,± = e- Q(T+)/2 e^ r+ 0_ for to > (23) 

These are 'normalized' to 2n5{u) — uj') with respect to (17), but are not yet defined in 
general over the whole spacetime. We wish to write the field operator ip(x) in terms of 
these modes as: 

^( x ) =ft+\ x ) +^ ( ~ ) (x) + ^ {0) (x) (24) 
where ^ {+ \x) = f j~ ^2 ^/W>+ 5 ^ 

J ~— TP D 



u=F,B 
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Figure 4. Representation of four of the nine possibilities for the structure of the causal 
envelope of an observer in 1+1 dimensional flat space. The dark curved lines represent 
a typical observer from each class, and the dotted lines represent typical 'hypersurfaces 
of simultaneity' for such observers. 



and ip(°\x) has compact support outside the causal envelope of our observer. To describe 
this construction, acceleration horizons must be considered. Consider the behavior of 
x ± (t\) as T\ — > oo. There are three possibilities: 



X 


'to - 




and 


x~ 


'(rx) ~ 


-> oo 


(A) 


x~ 


(rx) ~ 


-> X~ < OO 


and 


X 4 


(rx) ~ 


-> oo 


(B) 


X 4 


"(ta) - 


-> oo 


and 


x~ 


(rx) ~ 


-> oo 


(C) 



(x 4 ^ cannot both be finite, since ^— ^— = 1 for all t\). Similarly, as t\ — > — oo we have 
the possibilities 



x H 


-fa) - 


-> x 4 :^ > -oo 


and x~ 




-> — oo 


(A' 


x~ 


•(r A ) - 


-> xl^ > -oo 


and x 4 " 


(ta)- 


-> — oo 


(B' 


X H 




-> — oo 


and x~ 


(ta)- 


-> — oo 


(C 



These combine to give nine possible structures of the causal envelope of the observer. 
We consider four of these cases here (the other five are similar): 



• (C C) There are no horizons in this case. It is easy to see that the modes Vv°-,± 
are complete, so that we can achieve (24) with ?// )(x) = 0. 
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• (B C) In this case there is a future acceleration horizon at x~ = x^, which separates 
the spacetime into two regions - denoted I and II. Region I is the causal envelope of 
the observer, and region II the remainder. The backward-moving modes are defined 
naturally on the whole spacetime, since they are functions of x + only (the range of 
x + (t\) covers the whole real line). They form a complete set of backward moving 
modes. Although they were normalized on a hypersurface which is not Cauchy, it 
is straightforward to verify that the same norm holds on any Cauchy surface. The 
forward moving modes, on the other hand, are defined only in Region I. They can 
be extended to the whole spacetime by specifying i/) Ut F,±(x) — throughout region 
II. This ensures that the Dirac equation is satisfied even on the horizon (a wave- 
packet formed from these will be continuous on the horizon) and ensures that they 
are normalized to (2tc)5(uj — uo') on any Cauchy surface. Although they do not form 
a complete set, they can easily be supplemented by further 'forward moving modes', 
^i,ii{ x )i which are zero outside region II, and which together with the {i>u>,F,±(x)} 
comprise a complete set. (The ipi t n(x) might be chosen to be the forward moving 
modes of an observer in Case (C B'), which has x^ as a past horizon.) We can 
then achieve (24) with 

^°)( x) = J2^ n (x)c t 

i 

where the operators q satisfy anticommutation relations appropriate to the 
normalization of the There is considerable freedom in choosing these modes, 
which do not appear in the Hamiltonian H{r) (they all satisfy H\i^i : ii = 0) or 
in the number density operator N(x) (defined below) for this observer. These 
modes are 'unobservable' by our present observer, since they vanish throughout the 
corresponding causal envelope. 

• (B B') In this case there are future and past acceleration horizons at 

which divide the spacetime into regions I, II and III as shown in Figure 4. 
The backward-moving modes again form a complete set, defined on the whole 
spacetime. The forward- moving modes can be extended by specifying ipuj,F,±{z) = 
throughout regions I and III. They may be supplemented to form a complete set 
by including modes {ipij(x)} and {ipi>jn{x)}, with compact support in regions I 
and III respectively. Equation (24) is then achieved with 

ft°\x) = ^2^i,i(x)cij + ^2^i',iii(x)c i - JII 

i i' 

• (B A') This includes the well-studied uniform acceleration case (Section 2.2). It 
involves a future horizon at x~ = x^ and a past horizon at x + = xt^,, which 
intersect at a point, and divide the spacetime into 4 regions, denoted L, R, F and 
P in Figure 4. All hypersurfaces S r terminate at this point of intersection. In 
this case the backward-moving modes are defined in Regions R and F, and can be 
extended to the whole spacetime by specifying ipuj,B,±{x) = throughout regions 
L and P. They then satisfy (tp u ',B, \' I^Vb, a) = (^)S\yS(u — uj') on any Cauchy 
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surface. They can be supplemented by backward- moving modes ipi t Lp{x) having 
compact support in the union of Regions L and P. The forward-moving modes 
can be similarly extended, and supplemented by modes having compact 

support in the union of Regions L and F. Equation (24) is then achieved with 

ft°\ X ) = ^2^i, L p(x)Ci,LP + ^2^i>,LF(x)Ci>,LF 

i V 

Equation (24) has now been confirmed in all cases, and it is routine to substitute i\){x) 
into equation (15) to obtain the (pre-normal-ordered) Hamiltonian: 

^) = I E - (26) 

J a=F,B 

which has successfully been diagonalised, and does not depend on tp(°\x). The time- 
independence of the operators a Wi(7 , and of the Hamiltonian H{r) follows because 
the solutions ^ w ^±{x) are eigenstates of H\{t) for all r (and with the same eigenvalue 
oo for all r). 

The (second-quantized) total number operator N(t) [19, 25] can be written in terms 
of particle/antiparticle density operators A r± (x) as 

/oo 
dp{N + {x) + N~{x)) (27) 
-oo 

where N + (x) = ^ (+)t (x)M^ (+) (x) (28) 

oo oo 

= / I ^ W '^^- + ^-""*«Ub^b} (29) 



= 7V+(r-) + 7V+(r + ) (30) 
(say), and N~(x) = - : ^~ ) \x)M^-\x) : (31) 

oo oo 

= / £ I ^''"'""'^A- + ^ ,T %J^) (32) 



= N^(r-) + N B (r + ) (33) 

where the normal ordering is with respect to the observer's particle interpretation (i.e. 
the 6 Wi f)- By equating expressions (12) and (24) for ^(x), we find that 



—{{il) w ,o,+ \ u P,°,+) a P,« + (^,a,+ \u p , a ,-) bl a } (34) 

^{(^,a,-\Up,a,+)* al >a + (^ W)ff ,-|Up 1<T) -)*&p 1<y } ( 35 ) 

(A similar expression exists for the q.) Substitution into allows us to evaluate the 
particle content of the Dirac Vacuum \0m) as measured by this observer. For instance, 
the total number of particles/antiparticles measured is: 
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oo 

iv+ = (o M |iv + |o M ) = £ (0U,a 

oo 

TV- = (o M |iv-|o M > = IyY. (ry f W ( 36 ) 



<j=F,B 



where 



oo 

/W = <^,f,+ Kf,-> = f dpe- a ^ e ^P e ^-Ap) 

— oo 
oo 

&*,b = <^,b,+Kb,-> = J d P e a ^' 2 e^e^ x ti P ) (37) 

—00 

and 7^ = (^ CJ)0 ._|tt PiCT)+ ) = (3* pcT (38) 

Here /3 W p iCT and 7 wPi(7 are independent of r, since the states i^u>,a,± are all solutions of 
the Dirac equation. The variable p serves simply as an integration dummy. Though its 
physical origins come from an integral over space we could equally well use the label t\, 
viewing the integral as an integral along the observer's worldline. 

From equation (38) we see that (77^)^,0- = (PP^)u>u>,a for both a, so that from (36) 
the total number and frequency distribution of 'forward '-moving particles are equal to 
those of 'forward '-moving antiparticles, and similarly for 'backward'-moving particles. 
In general, the forward-moving particle content need not resemble the backward-moving 
particle content. 

The formula for the distribution of particles also takes a simple form: 

00 00 

n+(r,p) = (0 M \N+(x)\0 M ) = J^J e^"^ - 



00 00 

n F (r,p) ee (0 M \N-(x)\0 M ) = J^J e^'^" (39) 



= n%{r,p) (40) 

where equation (38) has been used to reach the last line. The spatial distributions of 
particles and antiparticles are the same, as expected from local conservation of charge. 
Further, the integral over p > which is implicit in the definition of (/3f3') uu i >er can be 
evaluated explicitly to yield: 

r\ POO pOO 

(A3 f W,F =- / dr a e- 2 ™ / dT d sm(2u a r d )g + (r a ,r d ) (41) 

J -00 JO 

where u a — + u/) uj d = — cu) (42) 

, , x 1 exp[^(a(r a + r d ) + a(r a - r d ))] 

and g±{T a ,T d ) = - 7^— ■ — r= (43) 

2r d J_ d dr exp [^a r a + r 
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In equation (39) this yields: 

±/ n 1 P dr , s f°°duj a 

n$(T ) = 4^Lt'^ T +T ' T) = L ^ jT ) (44) 

where 

2 f°° , sin[2cj(r- - r )l f 50 , , , , , ,„ r , 

«f )W (^ ) = — / dr a / dr d sm(2wr d )^ + (r , r d ) (45) 

ft J -oo T — T a Jq 

This density is a function only of t~ = r — p, as expected for forward- moving 
massless particles. It is defined such that n^(r~)dp gives the number of particles within 
dp of the point (t, p). Integration of n FjU (r) over p at any time r gives the total number 
of forward-moving particles of frequency u: 

/oo i*oo 
dpn FjU ,(T~) = / dT~n FjU (T~) (46) 
-oo J —oo 

In combination with equation (44), this result suggests that n F>(JJ {r~) should be 
interpreted as the frequency distribution of forward moving particles/antiparticles at 
the point r~ . Though this interpretation is reasonable, np tU1 (r~) must be averaged over 
a suitable range of uo and p before its predictions become more significant than the 
fluctuations [25]. This limitation is manifest for instance, in the fact that np yLU (T^) is 
not positive definite, whereas n^(r~) and ((3ft) WilJ ,F are positive definite. 

Equation (45), together with (43), specifies the distribution np^{ T ~) anywhere 
in the spacetime directly in terms of the observer's rapidity. Similarly, the spatial 
distribution of backward-moving particles/antiparticles is given by: 



oo oo 

n ± B (r,p) = J t J ^ m ^' B (47) 



-4_(t+ + t,t) = / — n B ,Ur + ) (48) 



4tt 2 r ° v ' ' L 2tt 



where 



. 2 f°° , sin[2cc;(r + - rj] Z" 00 , . /n . . . , AnS 
n B Ar + ) = — dr a — ^ / dr d sin{2ur d )g_{r a , r d ) (49) 

ft J -oo T — T a Jq 

Clearly the backward moving particle content (which exactly matches the backward 
moving antiparticle content) is a function only of r + = r + p. 

Although the distributions of 'forward'- and 'backward'-moving particles are not 
generally the same, we might expect them to be related when the observer's trajectory 
is time-symmetric. To study this possibility, consider an observer whose trajectory 
satisfies the time-reversal symmetry 

x\(-T\) = X\ (t\) 

tx(-rx) = -t x (r x ) (50) 

It follows that a(—T\) = — a(r A ), which gives n F ^{r') = 71b )W (— t') for all r'. On the 
hypersurface r = 0, for instance, where = ±p, this implies that the distribution of 
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forward moving particles exactly matches that of backward moving particles. However, 
on the observer's worldline, the distribution of forward moving particles is the time- 
reverse of the distribution for backward moving particles. 



4. Examples 

4-1. Inertial Observer 

For an inertial observer, a is constant, which gives g±(r a ,T d ) = 0. The particle content 
is therefore everywhere zero, as expected. 



4-2. Constant Acceleration 

For a uniformly accelerating observer we have 

9+(r a , r d ) = g-(T a ,r d ) = , ° r 

2r d 2 smh (ar d ) 

which is independent of r a . Hence the forward and backward moving particledistribu- 
tions are both uniform w.r.t. p for all r, and the frequency distribution is everywhere 
given by: 

1 f°° ( 1 a \ 

n F ,u = n B ,oj = - dr d sin(2^r d ) — — (51) 

tt Jo \ r d smh {aT d )J 



1 



2tylj 

1 + e~ 



(52) 



which is a thermal spectrum at temperature T = 2^^) as expected [8, 3, 26]. 

4-3. A 'Smooth Turnaround' Observer. 
Consider again the observer with acceleration 

a(r A ) = acosh~ 2 ( — 
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By substituting the corresponding rapidity a{j\) = ar c tanh(r/r c ) into equation (43), 
we immediately obtain the spatial distribution of forward or backward moving particles. 
At time r = these distributions are equal, as a result of the time-reversal invariance 
of the observer's trajectory. They are shown in Figure 5 as a function of ap for ar c = 1 
(bottom curve), 3,10,30,100 and oo (top line). As t c increases the particle density 
increases, and approaches the spatial uniformity found for the 'constant-acceleration' 
(t c — > oo) limit. 




0.2 0.4 0.6 0.8 1 0.1 072 O 0.4 0.5 



Figure 6(A). np^ip) for the smooth Figure 6(B). np^(p) for the smooth 

turnaround observer, as a function of uo/a for turnaround observer, as a function oiuj/a for 

p = = t, and ar c = 3, 10 and 30. r, ar c = 10 and ap = ±10. 

Figure 6 shows the frequency distribution np^(p) of forward- moving particles as a 
function of uj/a for r = 0. (This also equals the backward-moving distribution ns^(p) 
there.) In Figure 6 (A) we have p = and ar c = 3, 10, 30 and oo. We see clearly that the 
distribution approaches the thermal form as r c increases, with only the low frequencies 
differing significantly from the constant acceleration case. In Figure 6 (B), ar c = 10 and 
ap = +10 (top curve) and —10 (bottom curve). Since Uf,b depend only on r + and t~ 
respectively, Figure 6 (B) also represents the distribution of forward/backward moving 
particles on the observer's worldline, at t\ = 10/ a. At this time the observer sees more 
backward moving particles than forward moving particles, while the forward/backward 
moving distributions are the opposite at r = —10/ a. 

4-4- Acceleration at Late Times. 

Consider again the observer of Section 2.4, whose worldline is shown in Figure 3, in 
the (x, t) and (p, r) planes. Figure 7 shows the spatial distribution of forward and 
backward moving particles on this observer's trajectory. On the trajectory the observer 
consistently sees more forward than backward moving particles. Since np,B depend only 
on r + and r~ respectively we can again extrapolate from Figure 7 to the whole (p, r) 
plane. Well into region F (see Figure 3 (B)) the forward/backward distributions are 
approximately uniform, in agreement with the predictions of the constant acceleration 
case. In region L there are many forward moving particles but few backward movers. 
In region R there are many backward moving particles but few forward movers, while 
in region P, where the observer is almost inertial, few particles of any variety will be 
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Figure 7. np(T)/a (top) and ns(r)/a (bottom) for the 'late-time acceleration' 
observer as a function of ar for p = and ar c = 1 . 

observed. The distributions np^r) of Figure 7 can be decomposed into their frequency 
components by plotting n B ^{r) and n F ^{r) for various values of lu. This is done in 
Figure 8 for u/a = 0.05 (A) and 0.5 (B). 




Figure 8(A). np^r) and iib^{t) for the Figure 8(B). uf^(t) and iib^{t) for the 

'late-time acceleration' observer as a function 'late-time acceleration' observer as a function 



of ar for p = 0, w/a = 0.05 and ar c = 1. of ar for p = 0, lo / a = 0.5 and ar c = 1. 

These plots clearly show that np^{r) and ns^ir) are not positive definite, despite 
the positive definiteness of np^ij)- Also included in Figure 8 are the predictions 
obtained from the naive 'instantaneous thermal spectrum' 2tt/(1 + e 27ra V a ( T )) that would 
result from substituting a(r) directly into equation (52). The high frequency modes, 
which probe only short distance scales, closely trace these naive predictions, whereas 
the low frequency modes take more time to adjust from the 'inertial' to the 'uniform 
acceleration' predictions. 

Figure 9 shows n FlJJ {j) and n B ^{r) as a function of u/a, for r = — 30/a,0 and 
30 1 a. At late times the spectrum is close to thermal. For r = neither the forward 
(Figure 9 (B)) nor backward (Figure 9 (A)) moving distributions resemble a thermal 
spectrum, even at large frequencies, and ^.^^(t) is not positive definite. At much earlier 
times both distributions are very small, and neither is positive definite. The limit as 
uj — ► of these distributions is independent of r, reflecting the fact that the uj = mode 
is a global feature. However, this limit is different for forward and backward moving 
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1.2 



Figure 9(A). ub^{t) for the 'late-time 
acceleration' observer, as a function of 
uo/a 1 for p = 0, ar c = 1 and ar = 
—30,0 and 30. Also included are the 
'instantaneous thermal spectrum' predictions 

(1 +e 2™/a(r))-l ( dottcd ). 




Figure 9(B). n (t) for the 'late-time 
acceleration' observer, as a function of 
oj/a, for p = 0, ar c = 1 and ar = 
—30,0 and 30. Also included are the 
'instantaneous thermal spectrum' predictions 
(1 + e ^/a(r)yi (dotted). 



particles. 



5. Discussion 



The approach developed in [17, 19, 25] has been used to evaluate the particle content of 
the massless Dirac vacuum for an arbitrarily moving observer in 1+1 dimensions. This 
method uses Bondi's 'radar time', which provides a foliation of flat space that depends 
only on the observer's trajectory. It agrees with proper time on the trajectory, and is 
single valued in the causal envelope of the observer's worldline. 

The particle content can be described in terms of 'forward moving' particles (with 
particle density nj(r, p) dependent only on r — p), 'backward moving' particles (with 
particle density n~^(r, p) dependent only on r + p), and their corresponding antiparticles 
(described by Hp and ri^). The calculation of n| F (r, p) is largely independent of 
the possible presence of horizons, although when horizons exist an observer cannot 
assign particle/antiparticle densities to regions of spacetime invisible to this observer. 
This is entirely reasonable, though it differs from other approaches [15]. The particle 
content is clearly different for different observers, but for any given observer the total 
number of forward/backward moving particles is independent of r. This result is specific 
to massless particles in 1+1 dimensions, and is a direct consequence of conformal 
invariance. The spatial distribution of forward moving particles always matches that of 
forward moving antiparticles, and likewise for backward moving particles. In general, 
there is no connection between the observed numbers of forward moving particles and 
backward moving particles. However, when the observer's trajectory is time-reversal 
invariant, the distribution of forward and backward moving particles must match at 
r = 0. Their distributions at other times is then specified by the fact that n^(r, p) 
depends only on r — p and that n^(r, p) depends only on r + p. 

These general results have been illustrated by four examples. We have confirmed 
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that inertial observers do not observe particles in the inertial vacuum, and that uniformly 
accelerating observers measure the well-known Fermi-Dirac spectrum [3, 4, 5, 6, 7, 8, 10, 
25, 26]. This thermal spectrum is spatially uniform in 'radar distance' p and is the same 
for forward and backward moving particles, as expected. An observer with acceleration 
a(r\) = a cosh~ 2 (r / r c ) has also been studied. For large r c , this observer increasingly 
resembles the uniformly accelerating observer, although no horizons are present. For 
\p\, |r| << r c the particle content measured by this observer closely resembles the 
uniformly accelerating case, particularly for high frequency modes, which probe only 
short distance scales. This sheds light on a common misconception - that the Unruh 
effect is "due to the presence of a horizon" . As a result of the time-reversal invariance 
of this example, we find n^(0,p) = n^(0,p) for all p as generally predicted, although 
np{T,p) ^ n%(r,p) for r ^ 0. 

The further example treated here is an observer with a(r\) = 1+e ° /Te) ■ This observer 
accelerates uniformly at late times, but is inertial at early times. The observed particle 
content closely resembles the uniform acceleration result in the causal future of the 
acceleration, but approaches zero at early times (and away from the horizon). The 
distributions of forward and backward moving particles now differ significantly, with 
n^(r, 0) > n^(r, 0) for all r. This intriguing result suggests an observer-dependent 
violation of discrete symmetries, which calls for further investigation. 

Since massless fermions in 1+1 dimensions are conformally invariant, and 
consequently lack the scale dependence of the massive case (provided for instance by 
the Compton scale 1/m). The results obtained here will not generalize immediately 
to massive particles, or to 3+1 dimensions. Further work is needed to investigate 
those cases in more detail, and to explore the possible observer dependence of discrete 
symmetries. 

6. References 

[1] W. G. Unruh, Phys. Rev. D 14 (1976), 870 - 892. 

[2] P.C.W. Davics, J. Phys. A. 8(4) (1975), 609 - 616. 

[3] M. Soffel, B. Miiller and W. Grcincr, Phys. Rev. D 22 (1980), 1935 - 1937. 

[4] R. J. Hughes, Annals. Phys. 162 (1985), 1 - 30. 

[5] P. Candelas and D. Dcutsch, Proc. R. Soc. London Ser. A. 362 (1978), 251 - 262. 

[6] B. R. Iyer and A. Kumar, J. Phys. A. 13 (1980), 469 - 478. 

[7] S. Hacyan, Phys. Rev. D. 33(12) (1986), 3630 - 3633. 

[8] S. Takagi, Prog. Thco. Phys. Supplement No. 86, 1986. 

[9] N.D. Birrell and P.C.W. Davies, Quantum Fields in Curved Spacetime (Cambridge University 
Press, 1982). 

[10] M. Horibe, Prog. Theor. Phys. 61(2) (1979), 661 - 671. 

[11] L. Sriramkumar and T. Padmanabhan, Class. Quantum Grav. 13 (1996), 2061 - 2079. 

[12] N. Sanchez, Phys. Lett. B. 105(5) (1981), 375 - 380. 

[13] N. Sanchez, Phys. Lett. B. 87(3) (1979), 212 - 214. 

[14] N. Sanchez, Phys. Rev. D. 24(8) (1981), 2100 - 2110. 

[15] Zhu Jian-yang, Bao Aidong and Zhao Zheng, Intl. J. Theoret. Phys. 34(10) (1995), 2049 - 2059. 

[16] N. Sanchez and B. F. Whiting, Phys. Rev. D. 34(4), (1986), 1056 - 1071. 



The Fermionic Particle Density of Flat 1+1 Dimensional Spacetime seen by an Arbitrarily Moving Observer 

[17] C. E. Dolby, 'A State-Space Based Approach to Quantum Field Theory in Classical Background 

Fields', Thesis, Geometric Algebra Research Group, Cambridge University. 
[18] C. E. Dolby and S. F. Gull, Am. J. Phys. 69, (2001), 1257 - 1261. 
[19] C. E. Dolby and S. F. Gull, Annals. Phys. 293 (2001), 189 - 214. 
[20] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation (Freeman 1973). 
[21] M. Pauri and M. Vallisneri, Found. Phys. Lett. 13(5) (2000), 401 - 425. 
[22] H. Bondi, Assumption and Myth in Physical Theory (Cambridge University Press, 1967). 
[23] D. Bohm, The Special Theory of Relativity (W. A. Benjamin, 1965). 
[24] R. D'Inverno, Introducing Einstein's Relativity (Oxford University Press, 1992). 
[25] C. E. Dolby and S. F. Gull, 'Radar Time and a State-Space Based Approach to Quantum Field 

Theory in Gravitational and Electromagnetic Backgrounds', preprint number gr-qc/0207046. 
[26] W. Greiner, B. Miillcr and J. Rafelski, Quantum Electrodynamics of Strong Fields, with an 

Introduction to Modern Relativistic Quantum Mechanics, section 21 (Springer- Verlag, 1985). 



